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^ ■ Abstract 

Os ; 

Trilinear couplings of the neutral CP-even Higgs bosons in the Minimal Super- 
r^ [ symmetric Standard Model (MSSM) can be measured through the multiple produc- 

fr I tion of the lightest CP-even Higgs boson (h) at high-energy e~^e~ colliders. This 

^_ ■ includes the production of the heavier CP-even Higgs boson (H) via e~^e~ — > ZH, in 

association with the CP-odd Higgs boson (A) in e^e~ -^ AH, or via e~^e~ — > Ve^eH, 
with H subsequently decaying through H ^ hh. These processes can enable one to 
measure the trilinear Higgs couplings Xnhh and Xhhh, which can be used to theo- 
t~^ . retically reconstruct the Higgs potential. We delineate the regions of the MSSM 

r^ : parameter space m which these trhinear Higgs couplings could be measured. 

o 

a^ ; 1 Introduction 

^ ; 

Ph, Supersymmetry is at present the only known framework in which the Higgs sector [T of 

OhI the Standard Model (SM), so crucial for its internal consistency, is natural [2j. The mini- 

fy \ mal version of the Supersymmetric Standard Model (MSSM) contains two Higgs doublets 

(i^i, H2) with opposite hypercharges: Y{Hi) = —1, Y{H2) = +1, so as to generate masses 
for up- and down-type quarks (and leptons), and to cancel gauge anomalies. After spon- 
taneous symmetry breaking induced by the neutral components of Hi and H2 obtaining 
5^ i vacuum expectation values, {Hi) = vi, {H2) = V2, tan/3 = ^2/^1, the MSSM contains two 

neutral CP-even {h, H), one neutral CP-odd (A), and two charged (-ff^) Higgs bosons 
^ . Because of gauge invariance and supersymmetry, all the Higgs masses and the Higgs 
couplings in the MSSM can be described (at tree level) in terms of only two parameters, 
which are usually chosen to be tan/5 and ttia, the mass of the CP-odd Higgs boson. 

In particular, all the trilinear self-couplings of the physical Higgs particles can be pre- 
dicted theoretically (at the tree level) in terms of rriA and tan j3. Once a light Higgs boson 
is discovered, the measurement of these trilinear couplings can be used to reconstruct 
the Higgs potential of the MSSM. This will go a long way toward establishing the Higgs 
mechanism as the basic mechanism of spontaneous symmetry breaking in gauge theories. 
Although the measurement of all the Higgs couplings in the MSSM is a difficult task, 
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preliminary theoretical investigations by Plelin, Spira and Zerwas |3j, and by Djouadi, 
Haber and Zerwas (DHZ) j3], of the measurement of these couplings at the LHC and at 
a high-energy e~^e~ linear collider, respectively, are encouraging. 

We have considered in detail ^ the question of possible measurements of the trilinear 
Higgs couplings of the MSSM at a high-energy e~^e~ linear collider. We assume that such 
a facility will operate at an energy of 500 GeV with an integrated luminosity per year of 
-Cint = 500 fb~^ [6 . (This is a factor of 10 more than the earlier estimate.) In a later 
phase one may envisage an upgrade to an energy of 1.5 TeV. 

The trilinear Higgs couplings that are of interest are Xuhh, ^hhh, and Xhaa, involving 
both the CP-even and CP-odd Higgs bosons. The couplings Xnhh and Xhhh are rather 
small with respect to the corresponding trilinear coupling A^Jfh i^ ^^e SM (for a given 
mass of the lightest Higgs boson rrih), unless m/j is close to the upper value (decoupling 
limit). The coupling XhAA remains small for all parameters. 

Throughout, we include one-loop radiative corrections [Zj to the Higgs sector in the 
effective potential approximation. In particular, we take into account the parameters A 
and /i, the soft supersymmetry breaking trilinear parameter and the bilinear Higgs(ino) 
parameter in the superpotential, respectively, and as a consequence the left-right mixing 
in the squark sector, in our calculations. We thus include all the relevant parameters of 
the MSSM in our study jH]. Related work has recently been presented by Dubinin and 
Semenov jS|. 

For a given value of rrih, the values of these couplings significantly depend on the soft 
supersymmetry-breaking trilinear parameter A, as well as on /x, and thus on the resulting 
mixing in the squark sector. Since the trilinear couplings tend to be small, and depend 
on several parameters, their effects are somewhat difficult to estimate. 

The dominant source of multiple production of the Higgs (h) boson, is through Higgs- 
strahlung of H, and through production of H in association with the CP-odd Higgs boson. 
This source of multiple production can be used to extract the trilinear Higgs coupling 
Xuhh- The non-resonant fusion mechanism for multiple h production, e~^e~ -^ UgU^hh, 
involves two trilinear Higgs couplings, Xnhh and Xhhh, and is useful for extracting Xhhh- 

2 The Higgs Sector of the MSSM 

At the tree level, the Higgs sector of the MSSM is described by two parameters, which 
can be conveniently chosen as ttla and tan/3 [Tj. There are, however, substantial radiative 
corrections to the CP-even neutral Higgs masses and couplings jZHHj- They are, in general, 
positive, and they shift the mass of the lightest MSSM Higgs boson upwards. 

The Higgs mass falls rapidly at small values of tan/5. Since the LEP experiments 
are obtaining lower bounds on the mass of the lightest Higgs boson, they are beginning 
to rule out significant parts of the small-tan (3 parameter space, depending on the model 
assumptions. ALEPH finds a lower limit of ruh > 72.2 GeV, irrespective of tan/3, and a 
limit of ~ 88 GeV for 1 < tan/3 ^ 2 ^Hj- We take tan/? = 2 to be a representative value. 



3 Trilinear Higgs couplings 



In units of 5'm^/(2cos0w) = (vSGi?)^/^/?!!, the relevant tree-level trilinear Higgs cou- 
plings are given by 



X^Hhh = 2 sin 2a sin(/9 + a) — cos 2a cos(/5 -|- a) 
Khh = 3cos2asin(/5 + a), 
'^hAA = COS 2/5 sin (/5 + a), 



(3.1) 
(3.2) 
(3.3) 



with a the mixing angle in the CP-even Higgs sector, which can be calculated in terms 
of the parameters appearing in the CP-even Higgs mass matrix. The dominant one-loop 
radiative corrections are proportional to {mt/mw)'^ [TTj . 
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Figure 1: Trilinear Higgs couplings Xuhh, Xhhh and Xhaa as functions of rrih for two values 
of tan/3: (a) tan/3 = 2.0, (b) tan/? = 5.0. Each coupling is shown for three cases of 
the mixing parameters: no mixing (A = 0, ;U = 0, solid), mixing with A = 1 TeV and 
/i = — 1 TeV (dotted), as well as A = 1 TeV and /i = 1 TeV (dashed). 



The trilinear couplings depend significantly on ttia, and thus also on nih. This is shown 
in Fig. m where we compare Xuhh) Xhhh and XhAA for three different values of tan/3, and 
the SM quartic coupling A^'^ (which also includes one- loop radiative corrections [T^). 

At low values of rrih, the MSSM trilinear couplings are rather small. For some value 
of rrih the couplings Xnhh and Xhhh start to increase in magnitude, whereas XhAA remains 
small. The values of rrih at which they start becoming significant depend crucially on 
tan/?. For tan/? = 2 (Fig. [T^) this transition takes place around rrih ~ 90-100 GeV, 
whereas for tan/? = 5 the critical value of rrih increases to 100-110 (see Fig. ^p). In this 
region, the actual values of Xnhh and Xhhh (for a given value of rrih) change significantly if 



A becomes large and positive. A non-vanishing squark-niixing parameter A is thus quite 
important. Also, for special values of the parameters, the couplings may vanish P^. 

To sum up the behaviour of the trilinear couplings, we note that Xuhh and Xhhh are 
small for rrih ^ 100-120 GeV, depending on the value of tan/3. However, as rrih approaches 
its maximum value, which is reached rapidly as rriA becomes large, mA ^ 200 GeV, these 
trilinear couplings become large. Thus, as functions of m^, the trilinear couplings Xuhh 
and Xhhh are large for most of the parameter space. We also note that, for large values 
of tan/3, Xnhh tends to be relatively small, whereas Xhhh becomes large, if also m^ (or, 
equivalently, rrih) is large. 

4 Production mechanisms 

The different mechanisms for the multiple production of the MSSM Higgs bosons in e~^e~ 
collisions have been discussed by DHZ. The dominant mechanism for the production of 
multiple CP-even light Higgs bosons (h) is through the production of the heavy CP- 
even Higgs boson H, which then decays via H —>■ hh. The heavy Higgs boson H can 
be produced by if-strahlung, in association with A, and by the resonant WW fusion 
mechanism. These mechanisms for multiple production of h 

e+e- ^ ZH,AH \ rr , , /. in 

+ - - H- r ' H ^hh, 4.1 

are shown in Fig. [2l All the diagrams of Fig. |2l involve the trilinear coupling Xnhh- 
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Figure 2: Feynman diagrams for the resonant production of hh final states in e^e colli- 
sions. 

A background to ()4.ip comes from the production of the pseudoscalar A in association 
with h and its subsequent decay to hZ 

e+e" ^hA, A^ hZ, (4.2) 

leading to Zhh final states. A second mechanism for hh production is double Higgs- 
strahlung in the continuum with a Z boson in the final state, 

e+e" ^ Z* ^ Zhh. (4.3) 

We note that the non-resonant analogue of the Feynman diagram of Fig. |2|d involves, 
apart from the coupling Xnhk, the trilinear Higgs coupling Xhhh as well. 
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Figure 3: Feynman diagrams for the non- resonant WW fusion mechanism for the pro- 
duction of hh states in e~^e~ colhsions. 



Finally, there is a mechanism of multiple production of the lightest Higgs boson 
through non- resonant WW fusion in the continuum (see Section 7): 

e+e- -^ VeVeW*W* -^ u^uM- (4.4) 

It is important to note that all the diagrams of Fig. |21 involve the trilinear coupling 
Xnhh only. On the other hand, the non-resonant analogue of Fig. |2b, and Fig. Efc involve 
both the trilinear Higgs couplings Xnhh and Xhhh- 

5 Higgs-strahlung and Associated Production of H 

The dominant source for the production of multiple Higgs bosons (h) in e~^e~ collisions is 
through the production of the heavier CP-even Higgs boson H either via Higgs-strahlung 
or in association with A, followed, if kinematically allowed, by the cascade decay H —>■ hh. 
The cross sections for these processes can be found in [T^ ITK] . 

In Fig.|3]we plot the cross sections for the e~^e~ centre-of-mass energies ^/s = 500 GeV, 
as functions of the Higgs mass mn and for tan/5 = 2.0. For large values of the mass m^ 
of the pseudoscalar Higgs boson, all the Higgs bosons, except the lightest one (h), become 
heavy and decouple [16] from the rest of the spectrum. 

At values of tan/? that are not too large, the trilinear Hhh coupling Xnhh can be 
measured by the decay process H — >• hh, which has a width proportional to A^^^^^. How- 
ever, this is possible only if the decay is kinematically allowed, and the branching ratio 
is sizeable. In Fig. El we show the branching ratios (at tan/3 = 2) for the main decay 
modes of the heavy CP-even Higgs boson as a function of the H mass. Apart from the 
hh decay mode, the other important decay modes are H —>■ WW*, ZZ*. We note that the 
couplings of H to gauge bosons can be measured through the production cross sections 
for e~^e~ — > UgUeH; therefore the branching ratio BR[H -^ hh) can be used to measure 
the triple Higgs coupling Xnhh- 

For increasing values of tan/3, the Hhh coupling gradually gets weaker (see Fig. HJ, 
and hence the prospects for measuring Xnhh diminish. This is also indicated in Fig. El 
where we show the H branching ratios for tan/? = 5. 

There is actually a sizeable region in the m^-tan/3 plane where the decay H ^ hh is 
kinematically forbidden. This is shown in Fig. where we also display the regions where 
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Figure 4: Cross sections for the production of the heavy Higgs boson H in e^e~ colhsions, 
and for the background process in which Ah is produced. Sohd curves are for no mixing, 
y4 = 0, /i = 0. Dashed and dotted curves refer to mixing, as indicated. 



the H -^ hh branching ratio is in the range 0.1-0.9. Clearly, in the forbidden region, the 
^Hhh cannot be determined from resonant production. 



6 Double Higgs-strahlung and Triple h Production 

For small and moderate values of tan j3, the study of decays of the heavy CP-even Higgs 
boson H provides a means of determining the triple-Higgs coupling Xuhh- In order to 
extract the coupling Xhhh, other processes involving two-Higgs (h) final states must be 
considered. The Zhh final states, which can be produced in the non-resonant double 
Higgs-strahlung e'^e' — > Zhh , could provide one possible opportunity, since it involves 
the coupling Xhhh- These non-resonant processes have also been investigated [HE]. 

We show in Fig. U\ the Zhh cross section, with m = 1 TeV. The structure around 
rrih = 70 GeV (in the case of no mixing) is due to the vanishing and near- vanishing of the 
trilinear coupling. 

In the case of no mixing, there is a broad minimum from rrih — 78 to 90 GeV, followed 
by an enhancement around rrih ~ 90-100 GeV. This structure is due to the vanishing of 
the branching ratio for H -^ hh, which is kinematically forbidden in the region ruh — 78- 
90 GeV, see Fig. IHl (this coincides with the opening up of the channel H -^ WW), followed 
by an increase of the trilinear couplings. This particular structure depends considerably on 
the exact mass values ttih and rrih. Thus, it depends on details of the radiative corrections 
and on the mixing parameters A and /i. 
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Figure 5: Branching ratios for the decay modes of the CP-even heavy Higgs boson H^ for 
tan/3 = 2.0 and 5.0. Sohd curves are for no mixing, A = 0, /x = 0. Dashed and dotted 
curves refer to mixing, as indicated. 



7 Fusion Mechanism for Multiple-/^ Production 

A two-Higgs [hh) final state in e"*"e~ coUisions can also result from the WW fusion mech- 
anism, which can either be a resonant process as in ()4.H) . or a non-resonant one like ()4.4p . 
Since the neutral-current couplings are smaller than the charged-current ones, the cross 
section for the ZZ fusion mechanism in ()4.ip and ()4.4|) is an order of magnitude smaller 
than the WW fusion mechanism, and is here ignored. 

The WW fusion cross section for e+e~ -^ Hve^e can be written as [Hj (see also [Sj) 
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(7.1) 



This cross section is plotted in Fig. |3] for the centre-of-mass energy ^/s = 500 GeV, and 
for tan/? = 2.0, as a function of mn- The resonant fusion mechanism, which leads to 
[hh] + [missing energy] final states is competitive with the process e~^e~ — > HZ — > [hh] 
+ [missing energy], particularly at high energies. Since the dominant decay of h will be 
into bb pairs, the if-strahlung and the fusion mechanism will give rise to final states that 
will predominantly include four 6-quarks. On the other hand, the process e^e~ -^ AH 
will give rise to six 6-quarks in the final state, since the AH final state typically yields 
three-Higgs h[hh] final states. 

Besides the resonant WW fusion mechanism for the multiple production of h bosons. 
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Figure 6: The region in the m^-tan/3 plane where the decay H -^ hh is kinematically 
forbidden is indicated by a sohd hne contour. Also given are contours at which the 
branching ratio equals 0.1 (dotted), 0.5 (dashed) and 0.9 (dash-dotted, at the far left). 



there is also a non-resonant WW fusion mechanism: 

UpUphh, 



e e 



(7.2) 



through which the same final state of two h bosons can be produced. The cross section 
for this process, which arises through WW exchange as indicated in Fig. El can be written 
in the "effective WW approximation" as 



a(e e 



VeVehh) 



d.^*„..W. 



(7.3) 



where r = 4m|/s. Here, the cross section is written as a WW cross section, at invariant 
energy squared s = xs, folded with the WW "luminosity" |18j: 



dL(x) 
dx 



^f'^vf 



^^) - (l + -)log--2(l-x) 



(7.4) 



where v"^ + a? = 2. 

The WW cross section receives contributions from several amplitudes, according to 
the diagrams (a)-(d) in Fig. El We have evaluated these contributions 0. 

Our approach differs from that of DHZ in that we do not project out the longitudinal 
degrees of freedom of the intermediate W bosons. Instead, we follow the approach of 
Ref. [19], where transverse momenta are ignored everywhere except in the W propagators. 
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Figure 7: Cross section o"(e+e -^ Zhh) as a function of m/j. The dotted curve is the 
resonant production, the dashed curve gives the decouphng hmit. 



We show in Fig.|Hlthe WW fusion cross section, at -y/i = 1.5 TeV, as given by Eqs. ()7.1|) 
and ()7.3p . with m = 1 TeV. The structure is reminiscent of Fig. [71 and the reasons for 
this are same. Notice, however, that the scale is different. 



8 Sensitivity to Xnhh and Xhhh 

In Fig. El we have indicated in the mA-tan/3 plane the regions where Xnhh and Xhhh might 
be measurable for ^/s = 500 GeV. We identify regions according to the following criteria 

(i) Regions where Xuhh might become measurable are identified as those where (t{H) x 
BR(i7 — > hh) > 0.1 fb (solid), with the simultaneous requirement of 0.1 < BR(if -^ 
hh) < 0.9 [see Figs. ISHEj- In view of the recent, more optimistic, view on the 
luminosity that might become available, we also give the corresponding contours for 
0.05 fb (dashed) and 0.01 fb (dotted). 

(ii) Regions where Xhhh might become measurable are those where the continuum WW -^ 
hh cross section [Eq. ()7.3|) ] is larger than 0.1 fb (solid). Also included are contours 
at 0.05 (dashed) and 0.01 fb (dotted). 

Such regions are given for two cases of the mixing parameters A and /i, as indicated. We 
have excluded from the plots the region where rrih < 72.2 GeV, according to the LEP 
lower bound '\TU'^. This corresponds to low values of m^- 
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Figure 8: Cross section a{e^e — > Vf.Vf.hh) (via VrVT fusion) as a function of m^j. The 
dotted curve is the resonant production, the dashed curve gives the decouphng hmit. 



With an integrated luminosity of 500 fb~^, the contours at 0.1 fb correspond to 50 
events per year. This will of course be reduced by efficiencies, but should indicate the 
order of magnitude that can be reached. 

At y/s = 500 GeV, with a luminosity of 500 fb~^ per year, the trilinear coupling Xnhh 
is accessible in a considerable part of the m^-tan/? parameter space: at rriA of the order 
of 200-300 GeV and tan/9 up to the order of 5. With increasing luminosity, the region 
extends somewhat to higher values of m^. The "steep" edge around ttla — 200 GeV (where 
increased luminosity does not help) is determined by the vanishing of BR(iJ -^ hh) , see 
Fig. El 

The coupling Xhhh is accessible in a much larger part of this parameter space, but with 
a moderate luminosity, "large" values of tan/3 are accessible only if A is small. 

It should be stressed that the requirements discussed here are necessary, but not 
sufficient conditions for the trilinear couplings to be measurable. We also note that there 
might be sizable corrections to the WW approximation, and that it would be desirable 
to incorporate the dominant two-loop corrections to the trilinear couplings. 

9 Conclusions 

We have presented the results of a detailed investigation jH] of the possibility of measuring 
the MSSM trilinear couplings Xnhh and Xhhh at an e~^e~ collider. Where there is an 
overlap, we have confirmed the results of Ref. [4j. Our emphasis has been on taking into 
account all the parameters of the MSSM Higgs sector. We have studied the importance 
of mixing in the squark sector, as induced by the trilinear coupling A and the bilinear 
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Figure 9: Regions where trilinear couplings Xnhh and X^hh niight be measurable at ^/s = 
500 GeV. Inside contours labelled Xnhh, cr{H) x BR{H -^ hh) > 0.1 fb (solid), while 
0.1 < BR(iJ — > hh) < 0.9. Inside (to the right or below) contour labelled Xhhh, the 
continuum WW — > hh cross section exceeds 0.1 fb (solid). Analogous contours are given 
for 0.05 (dashed) and 0.01 fb (dotted). Two cases of mixing are considered, as indicated. 



coupling /i. 

At moderate energies {^/s = 500 GeV) the range in the ?72^-tan/5 plane that is ac- 
cessible for studying Xuhh changes quantitatively for non-zero values of the parameters A 
and fi. As far as the coupling Xhhh is concerned, however, there is a qualitative change 
from the case of no mixing in the squark sector. If A is large, then high luminosity is 
required to reach "high" values of tan/5. At higher energies {^/s = 1.5 TeV), the mix- 
ing parameters A and fi change the accessible region of the parameter space only in a 
quantitative manner. 
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